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6.6 Rayleigh-Darcy (or Horton-Rogers-Lapwood) in- 
stability in a porous layer 

6-6-Lapwood.tex 

Nield & Bejan, Chapter 6 Convection in Porous Media 
Related: Rayleigh-Bernard Problem (Chandrsekhar, Chapter II, Hydrodynamic and Hy- 
dromagnetic Stability) 

If a layer of viscous fluid is heated from below, instability can occur and leads to convec- 
tion cells important in meteorology. (Rayleigh-Benard Problem). 

If a saturated porous layer is heated from below, similar instabilily and convection can occur. 
This is of basic interest to geothermal convection and is relevant to the complex problem of 
heat transport due to the burial of nuclear waste in mountains or in a seabed. 

To give some visual ideas of what can happen in porous media, we shall borrow some 
photographic evidence for the mathematically similar Rayleigh-Benard problem of a pure 
fluid layer heated form below. 
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Figure 6.6.1: A saturated porus medium in geothermal gradient 



Static state 



u = 




Hence 




(6.6.1) 



Static equilibrium: 




hence 



Ps^Po-po9l^nz + (^(^j-2z^^ (6.6.2) 
Consider the perturbed state of smah disturbances: 

u = + u', T^Ts + r, P^ps+p' (6.6.3) 

then 

V ■ u' = (6.6.4) 



= -Vp' - ^u' + (igpoT'li (6.6.5) 

K 

dT' 



{pC)m^ + {pC) fu' ■ VTs = KrnV^r (6.6.6) 
6.6.1 Non-dimensionalization 

. Let Km, Km — Km/{pC)f be the conductivity and diflFusivity of the mixture, and k the 
permeabihty. Define 

{x,y,z) ^ h{x*,y*,z*), 1*, 

u' ^ ^u*, r = Are, p' ^ ^p* (6.6.7) 

h k 
Then, after omitting * for brevity, we get 

V ■ u = (6.6.8) 

= -Vp-u + Raek (6.6.9) 



dt 



(6.6.10) 



where 

Ra — ^^^^^ Rayleigh number in a porous medium (6.6.11) 

is the Rayleigh number (ratio of buoyancy force to diffusive resistence) of the porous medium. 
In Benard's problem, Rayleigh number is defined as 

Ra — — —, Rayleigh number in a pure fluid (6.6.12) 

fiK, 

Sometimes one calls 

k 

D = — , Darcy number (6.6.13) 
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Darcy number so that Rayleigh number of a porous medium is the product of the traditional 
Rayleigh number and Darcy number. 
Taking the curl of (6.6.9), 

V X u = Ra{i9y - j^^) (6.6.14) 

Note that the z compoent of the vorticity vector is zero. 
Taking the curl again and using 

V X V X u = V(V • u) - V^u 

we get 

Taking the z component, we get 



V^u; = i?aV^^ (6.6.15) 



where 

is the horizontal Laplacian. 

Equations (6.6.10) and (6.6.15) couple the two unknowns w and 9. The boundary con- 
ditions are 

w^e^O, z^O,l (6.6.17) 
After they are solved the other velocity components and pressure can be found. 

6.6.2 Solution for sinusoidal disturbances 

Let 

{w, 6) = {W{z), Q{z)) exp{ilx + imy - iujt) (6.6.18) 
and ^ 

dz 

then from (6.6.10), 

-iujQ -W ^{D'^ - a^)Q (6.6.19) 

and from (6.6.15) 

(D^ _ a^)w = -a^RaQ (6.6.20) 

where 

a^^e+m^ (6.6.21) 

The boundary conditions are : 

VF = e = 0, z = 0,l (6.6.22) 
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Because the equations and the boundary conditons are homogeneous, the problem for W 

and 6 is an eigenvalue problem. 

Note that I, m, a are related to dimensional wave numbers by 

I J 7 27r/i , , 27r/i , , „ , / 1 1 \ . , , fTZ TT 

I = kxh = ——, m = kyh = ——, a = kh ^ 2nh + with k = ^ kt. + kf, 

Ly \Ll LIJ y - y 

(6.6.23) 

6.6.3 Principle of exchange of stabilities 

We shall first show that oj must be purely imaginary. Multiplying (6.6.20) by W* and 
integrating from z = to z = 1, we get, after partial integration and using the boundary 
conditions, 

-f {\DW\'^ + a'^\W\'^)dz^-a^Ra f W*edz (6.6.24) 
Jo Jo 

Similarly we multipy (6.6.19) by G* and integrating from z = to z = 1, and get 

- / {\DQ\'^ + a-\e\^)dz ^ -ito [\Q\'^dz- [ We*dz (6.6.25) 
Jo J Jo 

Taking the complex conjugate of the second equation 

-/ i\De\^ + a''\e\'')dz^iuj* f\Q\^dz- f W*edz (6.6.26) 
Jo J Jo 

Eqs (6.6.24) and (6.6.26) can be combined by eliminating the cross product terms, , 

- / {\De\^ + a^\e\^)dz = -iLO* [ \e\'^dz - [ {\DW\'' + a'\W\'^)dz (6.6.27) 

Jo J ^ Jo 

Since all integrals above are real, —icv* — —Ui — iuj. must also be real. We conclude that 

Ur = 0, hence — iu — —uji (6.6.28) 

Marginal stability (the threshold of instability) occurs at lo^ = = 0. If uji > 0, the 
static state is unstable; coi = 0, marginally stable; if < 0, stable. A problem where the 
eigenfrequency is real so that marginal instabililty occurs when u; = is said to obey the 
principle of exchange of stabilities. 

6.6.4 Solution to eigenvalue problem 

Consider the situation at marginal stability : a; = 0, 

{D^ - a^)W = -a^Rae (6.6.29) 
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(6.6.30) 


Eliminating 6, we get 

(L>2 - a^fW = a^RaW 


(6.6.31) 


subject to 

W = ^, D^W = ^, z = 0,l 


(6.6.32) 


Expanding (6.6.31) 

D^W - 2a^D^W + a^W = a^RaW 


(6.6.33) 


Clearly D^W = on z = 0, 1. Differentiating (6.6.33) twice we see that = on z = 0, 1. 
Repeating the process we find 


D(2'«) = 0, m= 1,2,3,- •• , on z = 0, 1 


(6.6.34) 


Therefore the eigensolution must be 




W ~ sin j'ttz 


(6.6.35) 


To satisfy (6.6.31) it is necessary that 




[.•2„2 , 212 

U TT +aj ^ for j = 1,2, 3..., 


(6.6.36) 


which is the eigenvalue condition. For any j', Ra becomes unbounded for both - 
—> Qo and is curve concave upward in the plane of (abscissa) vs, Ra (ordinate) 
The lowest threshold occurs at j = 1, and 


and 


dRa 




i.e., 

2 2 
a — TT 


(6.6.37) 


or 

Rac = 477^ = 39.48 


(6.6.38) 



6.6.5 Possible convection patterns 

This is similar to Benard's problem which has been exhaustively studied theoretically and 
experimentally. There are many possibilities. Let us consider the lowest mode only with 

2-Dimensinal Rolls : {£ — Tr,m — 0) 
Take 

w; = cos TTX sin TTZ (6.6.39) 

then from mass conservation, 

= — sin TTX cos TTZ (6.6.40) 
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The dimensionless wavelength is = Ly = 2. Along lines x = 0, ±n, n = 1, 2, 3, 4, u = 
but to 7^ 0. Along x = 0, ±2m, w > 0, hence fluid rises vertically. Along x = ±2m — 1, 
< hence fluids sinks vertically. Along z = and l,w = 0. On the bottom (z = 0), n > 
while on the top (z = 1), if < x < 1, 2 < x < 3, 4 < x < 5, ... . The streamlines are shown 
in Figure 6.6.2. 




- i ^ / c I 2. 



Figure 6.6.2: Rolls in a period 



Rectangular cells: {£ = m — Tr/y/2). 



TTX Try , N 

w — cos cos sm TTZ (6.6.41) 

V 2 V 2 



From the z component of the vorticity equation (6.6.14), 

dv du 



and from continuity 

du ^ dv dw 
dx dy dz 

By cross differentiation, we get 

d'^u d'^u d'^w TT^ . TTX Try 

—r + — - = -TT^ = 1= Sm — COS COSTTZ 

and 

d'^v d'^v d^w TT^ TTX Try 

TT^ + TT^ = - „ „ = ^ cos —= sm ^ cos TTZ 

dx^ dy^ dzdx 2^2 a/2 a/2 



These are easily solved to give 



and 



1 . TTX Try 

■u = ;= sm — = cos — ;= COS Trz 

V2 V2 V2 



1 Trx . Try 

f = 1= COS — ;= sin —= cos Trz 

a/2 \/2 \/2 



The streamlines in a horizontal plane is shown in Figure 6.6.4. 
Hexagonal cells: See Chandrasekhar. 



FIgurm RectanguSar ^nd squara emits in a pura 

fiu'id hsatfid from bsiow. 
(AdRptBd from Chandrasakhar, Hydpodyitamfs and 
HydromagnaliG Stabitiiy). 



